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We analyze the logarithmic corrections due to ferromagnetic impurity ending bonds of open spin 
1 /2 antiferromagnetic chains, using the density matrix renormalization group technique. A universal 
finite size scaling ~ L ^ L for impurity contributions in the quasi-degenerate ground state energy 
is demonstrated for a zigzag spin 1/2 chain at the critical next nearest neighbor coupling and the 
standard Heisenberg spin 1/2 chain, in the long chain limit. Using an exact solution for the latter 
case it is argued that one can extract the impurity contributions to the entropy and specific heat 
from the scaling analysis. It is also shown that a pure spin 3/2 open Heisenberg chain belongs to 
the same universality class. 



PACS: 75.10.Jm, 75.40Mg 

The logarithmic corrections due to marginally irrel- 
evant operatorsa complicate greatly the comparison of 
experiments and numerical simulation results for finite 
size systems with analytical calculations. However, a 
careful analysis in some specific cases may yield useful 
information on the low energy excitation spectrum and 
relevant physical quantities. In this Report, we consider 
the finite size scaling for open spin 1/2 antiferromagnetic 
(AF) chains with ferromagnetic (FM) coupling at the 
ending bonds. This system exhibits the same behavior 
as a Kondo imparity coupled ferromagnetically to a Lut- 
tinger liquid.cT □ In these systems the Kondo screening is 
not complete, and the ground state is quasi-degenerate, 
i.e. the level spacing is vanishing faster than 1/L, where 
L is the system size. Recently, some exactly solvable 
models belonging to this class have been found, and the 
impurity entropy as well as specific heat hasi-ibeen ob- 
tained through thermodynamic Bethe ansatza a. We will 
show in this report the universal behavior of logarithmic 
corrections for this class of systems and, using the ex- 
act solution, we will argue that one can extract from the 
scaling analysis the impurity contributions to the entropy 
and specific heat when the exact solutions are not avail- 
able. We will also show that a pure open spin 3/2 AF 
chain (without additional FM ending bonds) belongs to 
the same universality class. 

The impurity effects in spin 1/2 Heisenberg AF chains 
have been discussed extensively in Ref. |7[ The logarith- 
mic corrections to scaling functions for spin 1/2 Hei&etu 
berg chains have been discussed in recent papers. Ertil 
Those corrections are due to bulk marginally irrele- 
vant operators, although their manifestations depend on 
boundary conditions^ Instead, we will carry out a de- 
tailed numerical analysis of the finite size scaling for the 
ground state near-degeneracy and low energy spectrum 
of open s — 1/2 spin chains due to FM ending bonds. 
To the best of our knowledge, this issue has not been 
addressed numerically up to now. We first consider the 
following Heisenberg chain with next nearest neighbor 



coupling, or equivalently, a zigzag chain: 

L-2 L-3 
H = Si ■ Sj+X + J2c Sj • S.; + 2 + Hi mp , 

i=2 i=2 
Hi m p = —Si • S 2 — Sl-i ■ Sl, (1) 

where L is the chain length and Sj is the s — 1/2 spin on 
site i. We draw the system in Fig. [lj. The nearest neigh- 
bor coupling is set to J — 1 and thejiext nearest neighbor 
coupling is set to the critical valueEa J 2 = J 2c = 0.2411. 
The two ending bonds are FM J' = —1, and the end- 
ing sites are impurity spins, also s = 1/2. There are no 
logarithmic corrections due to zero initial bulk marginal 
coupling at the critical value J 2 = J^Jiij'E-a Therefore, 
all logarithmic corrections in our calculations are coming 
entirely from the boundary effects. 
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FIG. 1. Open zigzag spin chain with coupling J and J2 for 
nearest and next nearest neighbors, respectively. Two ending 
impurity spins couple to the bulk by a FM bond J' < 0. 

A pure open zigzag chain without Hi mp has the same 
low energy spectrum as the pure open s = 1/2 chain 
when < J 2 < J 2c . There is a unique ground state 
with total spin 5 = 0, and one first excited state of spin 
5 = 1, with excitation energy scaled as irv/L for finite 
size chains,Q where v is the spin velocity. For the zigzag 
chain shown in Fig. [l] with FM ending bonds, the end- 
ing spins are not fully screened, and there is an RKKY 
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Ej 1 is not 



coupling between them scaled as Jrkky(L) = L1 ° gL in 
the large L limit. The two ending impurity spins form a 
singlet and a triplet with energy spacing Jrkky(L). We 
identify the following low energy states: 

(1) Two quasi-degenerate ground states: One is com- 
posed of the bulk 5 — state and singlet impurity 
state with energy Eq ; the other is formed by the bulk 
5 = state and the triplet impurity state with energy 

. We take Jrkky(L) = E? - E§ as the definition of 
Jrkky(L). 

(2) Four quasi-degenerate first excited states with excita- 
tion energy scaled as ttv/L. One is composed of the bulk 
5=1 state and singlet impurity state with energy E\. 
The other three are formed by the bulk first excited 5 = 1 
state and triplet impurity state with energies E^ 1 , E^ 1 , 
and Eq 1 , for the total spin S = 2, 1, and 0, respectively. 
Due to the bulk 5=1 excitation propagating in-between 
the ending spins, the energy difference Ef 1 
the same as Jrkky(L). 

We use. density matrix renormalization group (DMRG) 
methodEj to calculate low energy levels for the above 
Hamiltonian. By keeping m — 150 states, the truncation 
error is as small as 10~ 7 . We study even length chains 
only. The low-lying excitation energies are plotted vs 
1 / log L in Fig. g. In Fig. ||a we see the ground state is 
degenerate at the scale of the graph. The first four exci- 
tation energies scale to ttv/L. They correspond to E^ 1 , 
E(, Ef 1 , and Eq 1 , respectively, from bottom up. We note 
the ratio between the two energy spacings Eq 1 — E{ ! and 
E( ! — E^ 1 is approximately 1:2, which indicates these 
excitations can be identified as due to coupling of the bulk 
5=1 state and impurity triplet state. The next group 
of six energy levels scales to 2wv/L. They are composed 
of the. impurity singlet /triplet states and the two bulk 
levelsQ scaled as 2-kv/L of spin 5 = and 5 = 1. We 
have drawn guiding lines to separate these two groups of 
lowest excitations. In Fig. |^b, we magnify the scale of the 
energy spacing for the quasi-degenerate ground states to 
show the scaling E° — Eq = Jrkky(L) ~ L - Since 
there are no bulk logarithmic corrections involved at the 
critical coupling Ji c for zigzag chains, the logarithmic 
term appears solely due to the Kondo impurity effect. 
If the ending bound coupling J' = 0, the end impurity 
spins are decoupled from the bulk and the ground state 
is exactly degenerate. The impurity has nonzero entropy 
at zero temperature. When J' < 0, the ground state 
and low energy spectrum have an asymptotic degeneracy 
and the energy difference between these quasi-degenerate 



states scales as 



L log L ' 



This has been very clearly seen 



for the ground state. The zero temperature entropy will 
change, as we shall argue, depending on the coefficient a 
in the ground state energy scaling Jrkky{L) = L t ° L . 
We will see that this logarithmic scaling behavior due 
to boundary spins is universal for systems composed of 
Kondo impurities fcrromagnetically coupled to Luttinger 
liquids, even when there are also logarithmic corrections 
due to bulk marginally irrelevant operators. 
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2. Excitation energies times the chain length, 
Eq)L, vs 1/logL are plotted in figure (a) for the 
zigzag chain shown in Fig. ^. Scaling for the quasi-degenerate 
ground state, (Ei — Eq)L vs 1/ log L, is plotted in figure (b) 
along with a fitting line. 
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FIG. 3. Excitation energies times the chain length, 
(Lf — Eq)L, vs 1/logL are plotted in figure (a) for spin 
1/2 Heisenberg AF chain with FM ending bonds. Scaling for 
the quasi-degenerate ground state, (E± — Eq)L vs 1/logL, 
is plotted in figure (b) along with a fitting line. 

We consider now an open spin 1 /2 AF Heisenberg chain 
with impurity ending bonds described by the Hamilto- 



H 



L-2 



Si • s 2 



>L-1 



(2) 



The nearest neighbor coupling is set to J = 1 and the FM 
coupling J' for the ending spin is set to J' = — 1. For 
a pure spin 1/2 open chain without FM impurity bonds, 
the ground state energy scales asH E — e^L + ei — ^g- [1 + 
6/log 2 (L) + ...], where eo is the site energy, e± is the 
boundary energy, and v is the spin velocity for spin 1/2 
Heisenberg chain. The logarithmic correction appears 
here due to the bulk marginally irrelevant operator. We 
will demonstrate the ground state energy has one more 
in its finite size scaling due to FM Kondo 



term , . r 

L log L 

coupling: 



E? = e L + ei 



TTV 

24L 



LlogL 



+ 



(3) 



We calculate the energy levels by using DMRG method 
for even length chains. We keep m — 200 states and 
the truncation error is of the order 10~ 9 . The ground 
state is the same as for Hamiltonian (Q), with energy 
E§. The energy levels Ei are also labeled the same 
way. We plot the excitation energies (Ef — Eq)L vs 
1/logL in Fig. §a. In Fig.3b, the scaling E° - E° ~ 
, P 6 , is exhibited. The first excited states are four-fold 

L log L 

degenerate as we analyzed before for the zigzag chains. 
We have drawn guiding lines in Fig. [|to group these first 
excited states together. For the low energy spectrum, the 
excitations are again composed of the combined impurity 
and bulk spin states. The logarithmic scaling behavior of 
a standard Heisenberg chain due to FM impurity bonds 
is identical to that of zigzag chain, as the correction due 



to the bulk marginally irrelevant operator is of higher 
order. 
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th total spin S = 1 E° , obtained from the exact 
L } * L along with the fitting line. 



with 
nB is 



On the other hand, the mo£lel(g) has been solved ex- 
actly using the Bethe Ansatz,cl we can extract its ground 
state energy from the exact solution. We calculate, using 
the Bethe Ansatz equations of Ref. 0, the energy E° of 
5=1 state_for system size up to more than four thou- 
sands sitesJlj (We have not yet obtained the energy for 
another degenerate ground state with lower energy E® 
from the Bethe ansatz equations.) Following Eq.(||), we 
plot Ei — cqL + i^j- vs LX l gL in Fig. |^, where the site en- 
ergy eo = 1/4 — log 2, and the spin velocity v — ir/2. We 



obtain L7p 



e L 



24L 



0.787984- 



by the least 



L log L 

square fitting. Based on Eq.(|3]), we have e\ — 0.787984. 



0.66 
L log L 



from the 



We obtain also the scaling Jrkky (L) 
exact solution, to be compared with an almost identical 
fitting of the DMRG data (coefficient o = 0.6 instead of 
0.66 in the scaling formula Jrkky (L) — Ll ^ gL )- There- 
fore, we conclude that the quasi-degenerate ground state 
and low-lying excitations of Luttinger liquids with FM 
Kondo impurities can be indeed described by this loga- 
rithmic scaling function. n 

Moreover, using the exact solution,El the impurity en- 
tropy and specific heat can be expressed in terms of 
the coefficient a in the above finite size scaling. If two 
systems have the same ground state degeneracy and 
the same low-energy excitation spectrum, the thermo- 
dynamic properties at very low temperatures should also 
be identical to each other. Using this argument, we can 
estimate the impurity entropy and specific heat from the 
finite size scaling analysis of the energy spectrum even 
in the case when the exact solution is not available (e.g. 
the zigzag chain considered earlier). Using the same ar- 
gument we will discuss another related system, an open 
spin 3/2 AF chain without additional FM ending bonds. 
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FIG. 5. Excitation energies times chain length are plotted 
vs 1/logL for open s = 3/2 chains. The energy times chain 
length difference between the quasi-degenerate ground states 
is fit by 3.9/ log L. 

The spin 3/2 chain has been shown to have the same 
low energy physics as for spin 1/2 chainO When the 
chain is open, there are.-effective edge s' = 1/2-spins left 
at the ends of the chain.ES It was shown earlienlj that the 
energy spacing between the two degenerate ground states 
scales as Ll ^ gL - However, the scaling of the low energy 
spectrum is very difficult to study and the logarithmic 
corrections are big.EZI Following the way of dealing with 
such edge spins for integer spin chains,E3 it was proposed 
to treat them as impurity spins Jerromagnetically cou- 
pled to the bulk spin excitations .113 Unlike the previous 
cases, there are no additional FM impurity bonds at the 
ends of spin 3/2 chains. Now we calculate the scaling 
of the ground state energy spacing Jrkky(L) = L1 ° gL 
by DMRG. We keep m = 1000 states in DMRG and 
the truncation error is 10 -6 . We calculate only a few 
low energy levels and plot the excitation energies times 
length vs 1/ log L in Fig. 0. We obtain the coefficient 
a = 3.9 in scaling Jrkky(L) — L t ° g L ■ With the known 
spin velocity v = 3.87 for spin 3/2 chain,0 we obtain a 
nonzero entropy at zero temperature. Its exact value can 
be expressed in term of a analytically when exact solu- 
tion is available. (Non-zero impurity entropy is also pre- 
dicted in Takhatajan-Babujian spin-3/2 chain with spin- 
1/2 boundary impurities, where Bethe ansatz solution is 
available. □) We hope such a nonzero entropy at zero tem- 
perature can be measured in quasi- one-dimensional spin 
3/2 materials such as CsVCl 3 (Ref. |l]), AgCrP 2 S 6 (Ref. 
H), etc. 

In summary, we have carried out a detailed analysis 
of the low energy spectrum of FM Kondo impurity in 
Luttinger liquids, or equivalently, an open spin 1/2 AF 
Heisenberg chain. We have shown this class of models 
has a universal logarithmic quasi-dcgcncracy in the low 
energy states due to the RKKY interaction between the 
unscreened edge spins. We argue that nonzero entropy 



at zero temperature can be obtained for FM Kondo im- 
purity system by studying the finite size scaling of the 
ground state energies. 
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